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Preface 



In a remarkable development, I have constructed a new definition for a 
self-delimiting universal Turing machine (UTM) that is easy to program 
and runs very quickly. This provides a new foundation for algorithmic 
information theory (AIT), which is the theory of the size in bits of 
programs for self-delimiting UTM's. Previously, AIT had an abstract 
mathematical quality. Now it is possible to write down executable 
programs that embody the constructions in the proofs of theorems. So 
AIT goes from dealing with remote idealized mythical objects to being 
a theory about practical down-to-earth gadgets that one can actually 
play with and use. 

This new self-delimiting UTM is implemented via software written 
in C that is based on the related and now largely obsolete software that 
I presented in my report "Exhibiting randomness in arithmetic using 
Mathematica and C," IBM Research Report RC-18946, 94 pp., June 
1993. In its turn this report was a reworking of the software for my book 
Algorithmic Information Theory, Cambridge University Press, 1987. 

Using this new software, as well as the latest theoretical ideas, it is 
now possible to give a self-contained "hands on" mini-course presenting 
very concretely my latest proofs of my two fundamental information- 
theoretic incompleteness theorems. The first of these theorems states 
that an A'^-bit formal axiomatic system cannot enable one to exhibit 
any specific object with program-size complexity greater than N + c. 
The second of these theorems states that an A^-bit formal axiomatic 
system cannot enable one to determine more than N -\-c! scattered bits 
of the halting probability VL. 

Most people beheve that anything that is true is true for a reason. 
These theorems show that some things are true for no reason at all. 
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i.e., accidentally, or at random. 

The latest and I believe the deepest proofs of these two theorems 
were originally presented in my paper "Information-theoretic incom- 
pleteness," Applied Mathematics and Computation 52 (1992), pp. 83- 
101. This paper is reprinted in my book Information- Theoretic Incom- 
pleteness, World Scientific, 1992. 

As is shown in this course, the algorithms considered in the proofs 
of these two theorems are now easy to program and run, and by looking 
at the size in bits of these programs one can actually, for the first time, 
determine exact vahies for the constants c and c'. 

It is my intention to use this approach and software in an intensive 
short course on the limits of mathematics that I will give at the Univer- 
sity of Maine in Orono in the summer of 1994. I also intend to write up 
this course as a new book. This research report, however, is intended 
to make these important new ideas and software widely available in 
preliminary form as soon as possible. 

I wish to thank Prof. George Markowsky of the University of Maine 
at Orono for his stimulating invitation to give a summer workshop. I 
also thank Prof. John Casti of the Santa Fe Institute, Prof. Carlton 
Caves of the University of New Mexico, and Prof. Cristian Calude of 
the University of Auckland for stimulating discussions. 

I am grateful to IBM for its enthusiastic and unwavering support of 
my research for a quarter of a century, and to my current management 
chain at the IBM Research Division, Marty Hopkins, Eric Kronstadt, 
and Jeff Jaffe. Finally I thank the RISC project group, of which I am 
a member, for designing the marvelous IBM RISC System/6000 work- 
stations that I have used for all these calculations, and for providing 
me with the latest models of this spectacular computing equipment. 

All enquires, comments and suggestions regarding this software 
should be sent via e-mail to chaitin at watson. ibm. com. 
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Part I 
Explanation 



The New Idea 



First of all, note that we use here the toy LISP from my monograph 
Algorithmic Information Theory, Cambridge University Press, 1987. 
(Reprinted with minor revisions thrice, lastly in 1992.) 

Here is a quick summary of this toy LISP. Each LISP primitive 
function and variable is as before a single character, but they have 
been changed from the IBM APL2 character set to the standard ASCII 
character set. These primitive functions, all of which have a fixed num- 
ber of arguments, are now ' for QUOTE (1 argument), . for ATOM 
(1 argument), = for EQ (2 arguments), + for CAR (1 argument), - 
for CDR (1 argument), * for CONS (2 arguments), & for LAMBDA 
and DEFINE (2 arguments), : for LET-BE-IN (3 arguments), / for 
IF-THEN-ELSE (3 arguments), , for OUTPUT (1 argument), ! for 
EVAL (1 argument), and ? for depth-limited EVAL (had 2 arguments, 
now has 3). The meta- notation " indicates that an S-expression with 
explicit parentheses follows, not what is usually the case in my toy LISP, 
an M-expression, in which the parentheses for each primitive function 
arc implicit. Finally the empty list NIL is written (), and TRUE and 
FALSE are 1 and 0. 

The new idea that goes beyond what is presented in my Cam- 
bridge University Press monograph is this. We define our standard 
self-delimiting universal Turing machine as follows. Its program is in 
binary, and appears on a tape in the following form. First comes a 
toy LISP expression, written in ASCII with 7 bits per character. Note 
that this expression is self-delimiting because parentheses must bal- 
ance. The TM reads in this LISP expression, and then evaluates it. As 
it does this, two new primitive functions @ and % with no arguments 
may be used to read more from the TM tape. Both of these functions 
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explode if the tape is exhausted, kilhng the computation. (9 reads a 
single bit from the tape, and % reads in an entire LISP expression, in 
7-bit character chunks. 

This is the only way that information on the TM tape may be 
accessed, which forces it to be used in a self-delimiting fashion. This 
is because no algorithm can search for the end of the tape and then 
use the length of the tape as data in the computation. If an algorithm 
attempts to read a bit that is not on the tape, the algorithm aborts. 

How is information placed on the TM tape in the first place? Well, 
in the starting environment, the tape is empty and any attempt to read 
it will give an error message. To place information on the tape, a new 
argument has been added to the primitive function ? for depth-limited 
evaluation. 

Consider the three arguments a, /3 and 7 of ?. The meaning of the 
first argument, the depth limit a, has been changed slightly. If a is a 
non-null atom, then there is no depth limit. If a is the empty list, this 
means zero depth limit (no function calls or re-evaluations). And an 
A'^-element list a means depth limit N. The second argument /? of ? is 
as before the expression to be evaluated as long as the depth limit a 
is not exceeded. The new third argument 7 of ? is a list of bits to be 
used as the TM tape. 

The value 1/ returned by the primitive function ? is also changed. 
u is a. list. The first element of is ! if the evaluation of (3 aborted 
because an attempt was made to read a non-existent bit from the TM 
tape. The first element of z/ is ? if evaluation of /? aborted because the 
depth limit a was exceeded. ! and ? are the only possible error flags, 
because my toy LISP is designed with maximally permissive semantics. 
If the computation (3 terminated normally instead of aborting, the first 
element of u will be a list with only one element, which is the result 
produced by the computation (3. That's the first element of the list 1/ 
produced by the ? primitive function. 

The rest of the value z/ is a stack of all the arguments to the primitive 
function , that were encountered during the evaluation of (5. More 
precisely, if , was called N times during the evaluation of /3, then u 
will be a list of -I- 1 elements. The N arguments of , appear in u in 
inverse chronological order. Thus ? can not only be used to determine 
if a computation /3 reads too much tape or goes on too long (i.e., to 
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greater depth than a), but ? can also be used to capture all the output 
that /3 displayed as it went along, whether the computation /3 aborted 
or not. 

In summary, all that one has to do to simulate a self-delimiting 
universal Turing machine U (p) running on the binary program p is to 
write 

?0' !y.p 

This is an M-expression with parentheses omitted from primitive func- 
tions, because all primitive functions have a fixed number of arguments. 
With all parentheses supphed, it becomes the S-expression 

(?0('(!(%)))p) 

This says that one is to read a complete LISP S-expression from the TM 
tape p and then evaluate it without any time limit and using whatever 
is left on the tape p. 

Two more primitive functions have also been added, the 2-argument 
function ~ for APPEND, i.e., list concatenation, and the 1-argument 
function # for converting an expression into the list of the bits in its 
ASCII character string representation. These are used for construct- 
ing the bit strings that are then put on the TM tape using ?'s third 
argument 7. Note that the functions # and % could be programmed 
rather than included as built-in primitive functions, but it is extremely 
convenient and much much faster to provide them built in. 

Finally a new 1-argument identity function ~ with the side-effect of 
outputting its argument is provided for debugging. Output produced 
by ~ is invisible to the "official" ? and , output mechanism. ~ is 
needed because ?Q!/^7 suppresses all output 6 produced within its depth- 
controlled evaluation of (3. Instead ? stacks all output 6 from within 
(3 for inclusion in the final value u that ? returns, namely v = (atomic 
error flag or (value of (3) followed by the output 9). 
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Course Outline 



The course begins by explaining with examples my toy LISP. Possi- 
bly the theory of LISP program-size complexity is developed a little, 
following my recent papers "LISP program-size complexity," Applied 
Mathematics and Computation 49 (1992), pp. 79-93, "LISP program- 
size complexity II," Applied Mathematics and Computation 52 (1992), 
pp. 103-126, "LISP program-size complexity III," Applied Mathematics 
and Computation 52 (1992), pp. 127-139, "LISP program-size complex- 
ity IV," Applied Mathematics and Computation 52 (1992), pp. 141-147. 
These papers are reprinted in my book Information- Theoretic Incom- 
pleteness, World Scientific, 1992. 

LISP program-size complexity is extremely simple and concrete. In 
particular, it is easy to show that it is impossible to prove that a self- 
contained LISP expression is elegant, i.e., that no smaller expression 
has the same value. To prove that an A^-character LISP expression is 
elegant requires a formal axiomatic system that itself has LISP com- 
plexity N. Also, LISP program-size complexity is subadditive, because 
expressions are self-delimiting and can be concatenated, and also be- 
cause we are dealing with pure LISP and no side-effects get in the way. 
Moreover, the probability I^lisp that a LISP expression "halts" or has 
a value is well-defined, also because programs are self-delimiting. Fi- 
nally, it is easy to see that the LISP program-size complexity of the 
first bits of the LISP halting probability ^Ilisp grows linearly with 
A^. Therefore to be able to determine the first bits of I^lisp requires 
a formal axiomatic system whose LISP complexity also grows linearly 
with A^. 

It should then be pointed out that LISP programs have severe 
information-theoretic limitations because they do not encode informa- 
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tion very efficiently in 7-bit ASCII characters subject to LISP syntax 
constraints. Arbitrary binary programs are denser and much better, 
but they should at least be kept self-delimiting. 

So next we define our standard self-delimiting universal Turing ma- 
chine U{p) using 

?o' !y.p 

as explained in the previous chapter. 
Next we show that 

H{x,y) <H{x) + H{y) + c. 

Here H{- • ■) denotes the size in bits of the smallest program that makes 
our standard universal Turing machine compute ■ ■ ■. Thus this inequal- 
ity states that the information needed to compute the pair (x, y) is 
bounded by a constant c plus the sum of the information needed to 
compute X and the information needed to compute y. Consider 

*!%*!%() 

This is an M-expression with parentheses omitted from primitive func- 
tions. With all parentheses supplied, it becomes the S-expression 

(*(!(%))(*(!(%))())) 

In fact, c is just 7 times the size in characters of this LISP S-expression, 
which is exactly 20 characters. Thus c = 7 x 20 = 140 bits! See 
univ.lisp. Note that in standard LISP this would be something like 

(CONS (EVAL (READ-EXPRESSION)) 
(CONS (EVAL (READ-EXPRESSION)) 
NIL)) 

which is much more than 20 characters long. 

Looking at binary strings x with size in bits, we next show that 

H{x) < 2\x\ +c 

and 

H{x) < \x\ + H{\x\) + c'. 
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As before, the programs for doing this are exhibited and run. Also c 
and c' are determined. See univ.lisp. 

Next we show how to calculate the halting probability fl of our stan- 
dard self-delimiting universal Turing machine in the limit from below. 
The LISP program for doing this, omega, lisp, is now remarkably clear 
and fast, and much better than the one given in my Cambridge Uni- 
versity Press monograph. (See also omega2.1isp and omegaS. lisp.) 
Using this, we show that if Qat denotes the first bits of the fractional 
part of the base-two real number Q, then 

Hi^N) >N-c. 

Again this is done with a program that can actually be run and whose 
size gives us a value for c. See omega4.1isp. 

Next we turn to the self-delimiting program-size complexity H{X) 
for infinite r.e. sets X, which is not considered at all in my Cambridge 
University Press monograph. This is now defined to be the size in bits 
of the smallest LISP expression ^ that executes forever without halting 
and outputs the members of the r.e. set X using the LISP primitive 
, . , is an identity function with the side-effect of outputting the value 
of its argument. Note that this LISP expression ^ is allowed to read 
additional bits or expressions from the TM tape using the primitive 
functions ® and % if ^ so desires. But of course $, is charged for this; 
this adds to ^'s program size. 

It was in order to deal with such unending expressions ^ that the 
LISP primitive function for time-limited evaluation ? now captures all 
output from , within its second argument /3. 

To illustrate these new concepts, we show that 

H{X nY)< H{X) + H{Y) + c 

and 

H{X U y) < H{X) + H{Y) + c' 

for infinite r.e. sets X and Y . As before, we run examples and determine 
c and d . See setsO.lisp through sets4.1isp. 

Now consider a formal axiomatic system A of complexity A^, i.e., 
with a set of theorems Ta that considered as an r.e. set as above has 
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self-delimiting program-size complexity H(Ta). We show that A has the 
following limitations. First of all, we show directly that A cannot enable 
us to exhibit a specific S-expression s with self-delimiting complexity 
H{s) greater than N -\- c. See godel . lisp and godel2 . lisp. 

Secondly, using the lower bound of iV — c on H{Q^) established in 
omega4. lisp, we show that A cannot enable us to determine not only 
the first N + c' bits of O, but any N + c' bits of Q, even if they are scat- 
tered and we leave gaps. (See godel3 . lisp.) In my Cambridge Univer- 
sity Press monograph, this took a hundred pages to show, and involved 
the systematic development of a general theory using measure-theoretic 
arguments. Following "Information-theoretic incompleteness," Applied 
Mathematics and Computation 52 (1992), pp. 83-101, the proof is now 
a straight-forward Berry-paradox-like program-size argument. More- 
over we are using a deeper definition of H(A) = H(Ta) via infinite 
self-delimiting computations. 

And last but not least, the philosophical implications of all this 
should be discussed, especially the extent to which it tends to justify 
experimental mathematics. This would be along the lines of the discus- 
sion in my talk transcript "Randomness in arithmetic and the decline 
and fall of reductionism in pure mathematics," Bulletin of the European 
Association for Theoretical Computer Science, No. 50 (June 1993), pp. 
314-328. 

This concludes our "hand-on" mini-course on the information- 
theoretic limits of mathematics. 



Software User Guide 



All the software for this course is written in a toy version of LISP. 
*.lisp files are toy LISP code. One C program, lisp.c, is also pro- 
vided; it is the LISP interpreter. In these instructions we assume that 
this software is being run in a UNIX environment. 

To run the LISP interpreter, first compile it. This is done using 
the command cc -0 -olisp lisp.c. This interpreter is about 80 
megabytes big. If this is too large, reduce the #def ine SIZE before 
compiling it. 

To run the LISP interpreter interactively, that is, with input and 
output on the screen, enter 

lisp 

As each M-expression is read in, it is written out, then converted to an 
S-expression that is written out and evaluated.^] 

To run a LISP program X.lisp with output on the screen, enter 

lisp <X.lisp 

To run a LISP program X . lisp with output in a file rather than on the 
screen, enter 

lisp <X.lisp >X.run 



^The conversion from M-to S-expression mostly consists of making all implicit 
parentheses explicit. 
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[[[ 

First steps with my new construction for 
a self -delimiting universal Turing machine. 
We show that 

H(x,y) <= H(x) + H(y) + c 
and determine c . 

Consider a bit string x of length |x| . 
We also show that 

H(x) <= 2|x| + c 
and that 

H(x) <= |x| + H(|x|) + c 
and determine both these c's. 
]]] 

[first demo the new lisp primitive fimctions] 

" ' (1234567890) ' (abcdefghi) 

@ 

?0 '0 'O 
?0 '@ '(1) 
?0 '0 '(0) 
?0 '© '(x) 

?o '**@()**@() '(10) 

?0 '**,(§()**,(§()() '(10) 

?0 '**,(§()**,(§()**, 0() '(10) 

#'a 

#' (abcdef) 
#> (12(34)56) 
?0 '7. '(110 0001) 
?0 '% '(110 0010) 
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?0 '% '(110 0011) 
?0 '% '(110 0100) 
?0 "A '(110 0101) 
?0 >l '(110 0110) 
?0 '% ' (110 0111) 
?0 '7. #'a 
?0 '% '(010 100) 
?0 '7o '(010 1001) 

?0 '7, '(010 1000 oil 0001 oil 0010 010 1001) 
?0 '7. '(010 1000 oil 0001 oil 0010 ) 
?0 '7. '(010 1000 oil 0001 oil 001 ) 
?0 '7. '(010 100 ) 
?0 '7. #' (abcdef) 
?0 '7. #'(12(34)56) 
?0 '*7.*7.() #'a #'b 
?0 

':(f) :x@ :y@ /=xy *x(f) () (f) 
' (0011001101) 

#':(f) :x(§ :y@ /=xy *x(f) () (f) 

[ Here is the self-delimiting universal Turing machine ! ] 

[ (with slightly funny handling of out-of-tape condition) ] 

& (Up) ++?0' !7.p 

[Show that H(x) <= 2|x| + c] 

(U 

- ,#,':(f) :x@ :y(S /=xy *x(f) () (f) 
' (0011001101) 

) 

(U 

- ,#,':(f) :x@ :y(S /=xy *x(f) () (f) 
' (0011001100) 

) 

(U 

~ '*!7o*!7o() [The length of this bit string is the] 
[constant c in H(x,y) <= H(x)+H(y)+c . ] 

- #':(f) :x(§ :y@ /=xy *x(f) () (f) 

- '(0011001101) 

' #':(f) :x@ :y@ /=xy *x(f) () (f) 
' (1100110001) 

) 
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[Size of list in reverse decimal!] 
& (Se) /.e() (I[J(S-e)) 
[Increment reverse decimal] 
& (In) /.n'(l) :d+n :r-n 

/=dO*lr 

/=dl*2r 

/=d2*3r 

/=d3*4r 

/=d4*5r 

/=d5*6r 

/=d6*7r 

/=d7*8r 

/=d8*9r 
[/=d9]*0(Ir) 
[Reverse list] 
& (Re) /.e() "(R-e)*+e() 

[Convert to binary and display size in decimal] 

& (Me) :e [J#[Je :f ,(R[,](Se)) e 

(M'a) 

(M'O) 

& (Dk) /=0+k *l(D-k) /.-k () *0-k [D = decrement] 

,(D,(D,(D,(D,'(001))))) 

(U 

" ,#,' [Show that H(x) <= | x | + H( | x | ) + c] 

: (Re) /.e() '■(R-e)*+e() [R = reverse ] 

: (Dk) /=0+k *l(D-k) /.-k () *0-k [D = decrement] 
: (Lk) /.k *(§(L(Dk)) [L = loop ] 

(L(R!y.)) 

" #' ' (1000) 
'(0000 0001) 

) 
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omega, lisp 



[[[[ Omega in the limit from below! ]]]] 

[Look at the behavior of typical 7-bit programs] 
?0 '!//.' (010 1000) [Ipar] 

?0'\ (010 1001) [rpar] 
?0' ! (Oil 0001) [1] 

[All strings of length k / with same length as k ] 

& (Xk) /.k'(()) (Z(X-k)) 

[Append and 1 to each element of list] 

& (Zx) /.x() *"+x'(0) *"+x'(l) (Z-x) 

(Z'((a) (b))) 

(X'O) 

(X'(l)) 

(X'(ll)) 

(x'diD) 

[Size of list in reverse binary] 

& (Se) /.e() (I(S-e)) 

[Increment reverse binary] 

& (Ix) /.x'(l) /=+xO *l-x *0(I-x) 

(S'O) 

(S'(a)) 

(SHab)) 

(S'(abc)) 

(S'(abcd)) 

[Pad X to length of k on right and reverse] 
& (Rxk) /.X /.k() *0(Rx-k) " (R-x-k) *+x() 
(R'(l)'(ll)) 
(R'(Ol)'(llll)) 
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(R> (0001) ' (1111 1111)) 

[Set of programs in x that halt within time k] 
& (Hxk) /.x() /,=0.+?k' !%,+x *+x(H-xk) (H-xk) 
(H '((111) (111 1111) (000) (000 0000)) '0) 
(H , *#'a *#'('(xy)) *#':(X)(X)(X) () '(111)) 
[For LISP omega must separate read & exec . ] 
& (Gxk) /.x() [version of H for lisp omega ] 
: e ?0'7o,+x [read expression from prog +x] 
[If read finished, evaluate exp for time k ] 
[with empty tape, so @ and % will fail! ] 
: V /.+e e ?k++e()[run for time k, no tape] 
/,=0.+v *+x(G-xk) [program +x halted ] 
(G-xk) [program +x didn't halt ] 
(G '((111) (111 1111) (000) (000 0000)) '0) 
(G , *#'a *#'('(xy)) *#':(X)(X)(X) () '(111)) 
(H , *-#'@'(l) *-#'%#' (ab) '(111)) 
(G , *~#'@'(1) *-#'%#' (ab) '(111)) 
[Omega sub k!] 

& (Wk) *0*" . (R, (S, (H, (Xk)k))k) 

(WO) 

(W'(l)) 

(W'(ll)) 

(w'diD) 
[[[[ 

(w (111 D) 
(w'dii 11)) 
(w'dii 111)) 
]]]] 

(w'dii 1111)) 
[[[[ 

(w'dii 1111 D) 
]]]] 
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[[[[ Omega in the limit from below! ]]]] 

[All strings of length k / with same length as k ] 

& (Xk) /.k'(()) (Z(X-k)) 

[Append and 1 to each element of list] 

& (Zx) /.x() *"+x'(0) *"+x'(l) (Z-x) 

[Size of list in reverse binary] 

& (Se) /.e() (I(S-e)) 

[Increment reverse binary] 

& (Ix) /.x'(l) /=+xO *l-x *0(I-x) 

[Pad X to length of k on right and reverse] 

& (Rxk) /.X /.k() *0(Rx-k) ~ (R-x-k)*+x() 

[Set of programs in x that halt within time k] 

& (Hxk) /.x() /=0.+?k' ly.+x *+x(H-xk) (H-xk) 

[Omega sub k!] 

& (Wk) *0*" . (R, (S, (H, (Xk)k))k) 

(W'O) 

(W'(l)) 

(W'(ll)) 

(W'(lll)) 

(W'dii D) 
(w'Ciii 11)) 

(W'Clll 111)) 

(w'dii 1111)) 
(w'dii 1111 D) 



27 



28 The Limits of Mathematics (in C) 



omegaS.lisp 



[[[[ Omega in the limit from below! ]]]] 

[All strings of length k / with same length as k ] 

& (Xk) /.k'(()) (Z(X-k)) 

[Append and 1 to each element of list] 

& (Zx) /.x() *"+x'(0) *"+x'(l) (Z-x) 

[Size of list in reverse binary] 

& (Se) /.e() (I(S-e)) 

[Increment reverse binary] 

& (Ix) /.x'(l) /=+xO *l-x *0(I-x) 

[Pad X to length of k on right and reverse] 

& (Rxk) /.X /.k() *0(Rx-k) ~ (R-x-k)*+x() 

[Set of programs in x that halt within time k] 

& (Hxk) /.x() /=0.+?k' ly.+x *+x(H-xk) (H-xk) 

[Omega sub k!] 

& (Wk) *0*" . (R(S(H(Xk)k))k) 

(W'O) 

(W'(l)) 

(W'(ll)) 

(W'(lll)) 

(W'dii D) 
(w'Ciii 11)) 

(W'Clll 111)) 

(w'dii 1111)) 
(w (111 1111 D) 

(W (111 1111 11)) 
(W (111 1111 111)) 

(w'dii 1111 111 D) 
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(111 1111 111 11)) 

(W'Clll 1111 111 111)) 

(w'dii 1111 111 1111)) 

(w (111 1111 111 1111 D) 

(W'Ciii 1111 111 1111 ID) 

[The following exhaust storage:] 

[[[[ 

(W (111 1111 111 1111 111)) 

(w'dii 1111 111 1111 111 D) 

(w'dii 1111 111 1111 111 ID) 

(w'dii 1111 111 1111 111 111)) 

(w'dii 1111 111 1111 111 iiiD) 
]]]] 



omega4.1isp 



[[[ 

Show that H(Omega sub n) > n - c and determine c. 
Omega sub n is the first n bits of Omega. 
]]] 

[First test new stuff] 

[Compare two fractional binary numbers, i.e., is 0.x < O.y ?] 
& (<xy) /.X /.y 

«'(0)y) 
/•y 

/+x /+y (<-x-y) 


/+y 1 

(<-x-y) 

«' (000) '(000)) 

(<' (000) ' (001)) 
(<' (001) ' (000)) 
(<' (001) ' (001)) 
«' (110) '(110)) 
(<'(110)'(111)) 
(<' (111) ' (110)) 

(<' (111) ' (111)) 

«'()'(000)) 

(<'()' (001)) 
(<'(000)'()) 
«'(001)'()) 
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[Now run it all!] 
++?0' !% 
> * > 

[All strings of length k / with same length as k ] 

: (Xk) /.k'(()) (Z(X-k)) 

[Append and 1 to each element of list] 

: (Zx) /.x() *"+x'(0) *"+x'(l) (Z-x) 

[Size of list in reverse binary] 

: (Se) /.e() (I(S-e)) 

[Increment reverse binary] 

: (Ix) /.x'(l) /=+xO *l-x *0(I-x) 

[Pad X to length of k on right and reverse] 

: (Rxk) /.X /.k() *0(Rx-k) " (R-x-k)*+x() 

[Set of programs in x that halt within time k] 

: (Hxk) /.x() /~=0.+?k' !7.~+x *+x(H-xk) (H-xk) 

[Omega sub k without 0. at beginning 

(i.e. only fractional part).] 
: (Wk) (R(S(H(Xk)k))k) 

[Compare two fractional binary numbers, i.e., is 0.x < O.y ?] 
: (<xy) /.X /.y 

«'(0)y) 

/•y 

/+x /+y (<-x-y) 


/+y 1 

«-x-y) 

: w !°/o [Read and execute from remainder of tape 

a program to compute an n-bit 
initial piece of Omega.] 
: (Lk) [Main Loop] 

: X (Wk) [Compute the kth lower bound on Omega] 

/(<xw) (L*lk) [Are the first n bits OK? If not, bump k.] 
(B(Xk)) [Form the union of all output of k-bit 
programs within time k, output it, 
and halt . 



omegaA.lisp 
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This is bigger than anything of complexity 
less than or equal to n!] 
[This total output will be bigger than each individual output, 
and therefore must come from a program with more than n bits . 

] 

[Union of all output of programs in list p within time k.] 
: (Bp) /.p() * ~?k'!yo~+p (B-p) [ k is implicit argument.] 

(L()) [Start main loop running with k initially zero.] 

'(1111) [These really are the first 4 bits of Omega!] 
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setsO.lisp 



[[[ 

Test basic (finite) set functions. 
]]] 

[Set membership predicate; is e in set s?] 
& (Mes) /.sO /=e+sl (Me-s) 
(Mx' (12345xabcdef)) 
(Mq' (12345xabcdef)) 

[Eliminate duplicate elements from set s] 
& (Ds) /.s() /(M+s-s) (D-s) *+s(D-s) 
(D' (1234512345abcdef)) 
(D(D' (1234512345abcdef))) 
[Set union] 

& (Uxy) /.xy /(M+xy) (U-xy) *+x(U-xy) 
(U' (12345abcdef) ' (abcdef UVWXYZ) ) 
[Set intersection] 

& (Ixy) /.x() /(M+xy) *+x(I-xy) (I-xy) 
C (12345abcdef) ' (abcdef UVWXYZ) ) 
[Subtract set y from set x] 
& (Sxy) /.x() /(M+xy) (S-xy) *+x(S-xy) 
(S' (12345abcdef ) ' (abcdef UVWXYZ) ) 

[Identity function that outputs a set of elements] 
& (Os) /.s() *,+s(0-s) 
(0' (12345abcdef)) 

[Combine two bit strings by interleaving them] 

& (Cxy) /.xy /.yx *+x*+y (C-x-y) 

(C (0000000000) ' (11111111111111111111)) 
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setsl.lisp 



[[[ 

Show that 

H(X set-union Y) <= H(X) + H(Y) + c 
and that 

H(X set-intersection Y) <= H(X) + H(Y) + c 
and determine both c's. 
Here X and Y are INFINITE sets. 
]]] 

[Combine two bit strings by interleaving them] 
& (Cxy) /.xy /.yx *+x*+y (C-x-y) 

[[[++]]]?0'!% 

',#,' 

[Package of set functions from setsO.lisp] 

: (Mes) /.sO /=e+sl (Me-s) 

: (Ds) /.s() /(M+s-s) (D-s) *+s(D-s) 

: (Uxy) /.xy /(M+xy) (U-xy) *+x(U-xy) 

: (Ixy) /.x() /(M+xy) *+x(I-xy) (I-xy) 

: (Sxy) /.x() /(M+xy) (S-xy) *+x(S-xy) 

: (Os) /.s() *~,+s(0-s) [<===cheating to get display!] 

[Main Loop: 

o is set of elements output so far. 
For first set, 

t is depth limit (time), b is bits read so far. 
For second set, 

u is depth limit (time), c is bits read so far. 
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] 

: (Lotbuc) 

V ~?~t'!yo~b [Run 1st computation again.] 
w ~?~u'!7o~c [Run 2nd computation again.] 
X (U-v-w) [Form UNION of sets so far] 

y (O(Sxo)) [Output all new elements] 

[This is an infinite loop. But to make debugging easier, 
halt if both computations have halted.] 
/ /=0.+v /=0.+w 100 X [If halts, value is output so far] 
[Bump everything before branching to head of loop] 
(L X [Value of y is discarded, x is new o] 

/="?+v *lt t [Increment time for 1st computation] 

/="!+v ~'b*®0 b [Increment tape for 1st computation] 
/="?+w *lu u [Increment time for 2nd computation] 

/="!+w "c*@() c [Increment tape for 2nd computation] 

) 

(L() ) [Initially all 5 induction variables 

are empty.] 



(C [Combine their bits in order needed] 

[Wrong order if programs use @ or %] 
,#,'*, a*, b*,cO [Program to enumerate 1st FINITE set] 

,#, '*,b*,c*,dO [Program to enumerate 2nd FINITE set] 

) 



sets2.1isp 



[[[ 

Show that 

H(X set-union Y) <= H(X) + H(Y) + c 
and that 

H(X set-intersection Y) <= H(X) + H(Y) + c 
and determine both c's. 
Here X and Y are INFINITE sets. 
]]] 

[Combine two bit strings by interleaving them] 
& (Cxy) /.xy /.yx *+x*+y (C-x-y) 

[[[++]]]?0'!% 

',#,' 

[Package of set functions from setsO.lisp] 

: (Mes) /.sO /=e+sl (Me-s) 

: (Ds) /.s() /(M+s-s) (D-s) *+s(D-s) 

: (Uxy) /.xy /(M+xy) (U-xy) *+x(U-xy) 

: (Ixy) /.x() /(M+xy) *+x(I-xy) (I-xy) 

: (Sxy) /.x() /(M+xy) (S-xy) *+x(S-xy) 

: (Os) /.s() *~,+s(0-s) [<===cheating to get display!] 

[Main Loop: 

o is set of elements output so far. 
For first set, 

t is depth limit (time), b is bits read so far. 
For second set, 

u is depth limit (time), c is bits read so far. 



39 



40 



The Limits of Mathematics (in C) 



] 

: (Lotbuc) 

V ~?~t'!yo~b [Run 1st computation again.] 
w ~?~u'!7o~c [Run 2nd computation again.] 
X (I-v-w) [Form INTERSECTION of sets so far] 

y (O(Sxo)) [Output all new elements] 

[This is an infinite loop. But to make debugging easier, 
halt if both computations have halted.] 
/ /=0.+v /=0.+w 100 X [If halts, value is output so far] 
[Bump everything before branching to head of loop] 
(L X [Value of y is discarded, x is new o] 

/="?+v *lt t [Increment time for 1st computation] 

/="!+v ~'b*®0 b [Increment tape for 1st computation] 
/="?+w *lu u [Increment time for 2nd computation] 

/="!+w "c*@() c [Increment tape for 2nd computation] 

) 

(L() ) [Initially all 5 induction variables 

are empty.] 



(C [Combine their bits in order needed] 

[Wrong order if programs use @ or %] 
,#,'*, a*, b*,cO [Program to enumerate 1st FINITE set] 

,#, '*,b*,c*,dO [Program to enumerate 2nd FINITE set] 

) 



setsS.lisp 



[[[ 

Show that 

H(X set-union Y) <= H(X) + H(Y) + c 
and that 

H(X set-intersection Y) <= H(X) + H(Y) + c 
and determine both c's. 
Here X and Y are INFINITE sets. 
]]] 

[Combine two bit strings by interleaving them] 
& (Cxy) /.xy /.yx *+x*+y (C-x-y) 

[IMPORTANT: This test case never halts] [<=====!!!] 

[[[++]]]?0'!% 

',#,' 

[Package of set functions from setsO.lisp] 

: (Mes) /.sO /=e+sl (Me-s) 

: (Ds) /.s() /(M+s-s) (D-s) *+s(D-s) 

: (Uxy) /.xy /(M+xy) (U-xy) *+x(U-xy) 

: (Ixy) /.x() /(M+xy) *+x(I-xy) (I-xy) 

: (Sxy) /.x() /(M+xy) (S-xy) *+x(S-xy) 

: (Os) /.s() *~,+s(0-s) [<===cheating to get display!] 

[Main Loop: 

o is set of elements output so far. 
For first set, 

t is depth limit (time), b is bits read so far. 
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For second set, 

u is depth limit (time), c is bits read so far. 

] 

: (Lotbuc) 

V Tt'l'/ob [Run 1st computation again.] 

w Tu'l'/oC [Run 2nd computation again.] 

X (I-v-w) [Form INTERSECTION of sets so far] 

y (O(Sxo)) [Output all new elements] 
[This is an infinite loop. But to make debugging easier, 
halt if both computations have halted.] 
/ /=0.+v /=0.+w 100 X [If halts, value is output so far] 
[Bump everything before branching to head of loop] 
(L X [Value of y is discarded, x is new o] 

/="?+v *lt t [Increment time for 1st computation] 

/="!+v "b*@() b [Increment tape for 1st computation] 
/="?+w *lu u [Increment time for 2nd computation] 

/="!+w "c*(9() c [Increment tape for 2nd computation] 

) 

(L() ) [Initially all 5 induction variables 

are empty.] 



(C [Combine their bits in order needed] 

[Wrong order if programs use or %] 
[Program to enumerate 1st INFINITE set] 

,#,':(Lk) (L,*l*lk) (LO) 

[Program to enumerate 2nd INFINITE set] 

,#,':(Lk) (L,*l*l*lk) (LO) 

) 



sets4.1isp 



[[[ 

Show that 

H(X set-union Y) <= H(X) + H(Y) + c 
and that 

H(X set-intersection Y) <= H(X) + H(Y) + c 
and determine both c's. 
Here X and Y are INFINITE sets. 
]]] 

[Combine two bit strings by interleaving them] 
& (Cxy) /.xy /.yx *+x*+y (C-x-y) 

[IMPORTANT: This test case never halts] [<=====!!!] 

[[[++]]]?0'!% 

',#,' 

[Package of set functions from setsO.lisp] 

: (Mes) /.sO /=e+sl (Me-s) 

: (Ds) /.s() /(M+s-s) (D-s) *+s(D-s) 

: (Uxy) /.xy /(M+xy) (U-xy) *+x(U-xy) 

: (Ixy) /.x() /(M+xy) *+x(I-xy) (I-xy) 

: (Sxy) /.x() /(M+xy) (S-xy) *+x(S-xy) 

: (Os) /.s() *~,+s(0-s) [<===cheating to get display!] 

[Main Loop: 

o is set of elements output so far. 
For first set, 

t is depth limit (time), b is bits read so far. 
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For second set, 

u is depth limit (time), c is bits read so far. 

] 

: (Lotbuc) 

V Tt'l'/ob [Run 1st computation again.] 

w Tu'l'/oC [Run 2nd computation again.] 

X (U-v-w) [Form UNION of sets so far] 

y (O(Sxo)) [Output all new elements] 
[This is an infinite loop. But to make debugging easier, 
halt if both computations have halted.] 
/ /=0.+v /=0.+w 100 X [If halts, value is output so far] 
[Bump everything before branching to head of loop] 
(L X [Value of y is discarded, x is new o] 

/="?+v *lt t [Increment time for 1st computation] 

/="!+v "b*@() b [Increment tape for 1st computation] 
/="?+w *lu u [Increment time for 2nd computation] 

/="!+w "c*(9() c [Increment tape for 2nd computation] 

) 

(L() ) [Initially all 5 induction variables 

are empty.] 



(C [Combine their bits in order needed] 

[Wrong order if programs use or %] 
[Program to enumerate 1st INFINITE set] 

,#,':(Lk) (L,*l*lk) (LO) 

[Program to enumerate 2nd INFINITE set] 

,#, ' : (Lk) (L,*2*2*2k) (L()) 

) 



godel.lisp 



[[[ 

Show that a formal system of complexity N 
can't prove that a specific object has 
complexity > N + c, and also determine c. 
Formal system is a never halting lisp expression 
that output pairs (lisp object, lower bound 
on its complexity). E.g., (x(llll)) means 
that X has complexity H(x) greater than 4. 
]]] 

[ (<xy) tells if x is less than y ] 
& (<xy) /.X /.yOl 

/.yO«-x-y) 

«'(111)'(11)) 

[ Examine pairs in p to see if 2nd is greater than n ] 

[ returns to indicate not found, or pair if found ] 

& (Epn) /.p /(<n+-+p) +p (E-pn) 

(E'((x(ll))(y(lll)))'()) 

(E'((x(ll))(y(lll)))'(l)) 

(E'((x(ll))(y(lll)))'(ll)) 

(E'((x(ll))(y(lll)))Hlll)) 

(E'((x(ll))(y(lll)))'(llll)) 

++?0 ' ! % 

> ''^ > 
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[ (<xy) tells if x is less than y ] 
: (<xy) /.X /.yOl 

/.yO«-x-y) 
[ Over-write real definition for test ] 
: (<xy) 1 

[ Examine pairs in p to see if 2nd is greater than n ] 
[ returns to indicate not found, or pair if found ] 
: (Epn) /.p /(<n+-+p) +p (E-pn) 

[Parameter in proof] 



k 


(11111) 


k 


~ 'kk 


k 


~ "kk 


k 


~ ~kk 


k 


~ ~kk 



[Main Loop - t is depth limit (time) , b is bits read so far] 
: (Ltb) 

: V ~?~tM°/o~b [run universal computer again] 

: s (E-v"kb) [look for pair with 2nd > 16k + # of bits read] 

/s +s [Found it ! Output 1st and halt] 

/="!+v (Lt"b*@()) [Read another bit from program tape] 

/="?+v (L*ltb) [Increase depth/time limit] 

"? [Surprise, formal system halts, so we do too] 

(L()()) [Initially, depth limit and no bits read] 

[ 

,#,','((xy)(ll)) 
] 

,#,','(x()) 



godel2.1isp 



[[[ 

Show that a formal system of complexity N 
can't prove that a specific object has 
complexity > N + c, and also determine c. 
Formal system is a never halting lisp expression 
that output pairs (lisp object, lower bound 
on its complexity). E.g., (x(llll)) means 
that X has complexity H(x) greater than 4. 
]]] 



++?o ' ! y. 



>^ > 

[ (<xy) tells if x is less than y ] 
: (<xy) /.X /.yOl 

/.yO(<-x-y) 



[ Examine pairs in p to see if 2nd is greater than n ] 
[ returns to indicate not found, or pair if found ] 
: (Epn) /.p /(<n+-+p) +p (E-pn) 



[Parameter in proof] 



k ~ 


(11111) 


k ~ 


~kk 


k ~ 


~kk 


k ~ 


~kk 


k ~ 


'kk 



[Main Loop - t is depth limit (time) , b is bits read so far] 
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: (Ltb) 

: V ~?~t' \y,~'b [run universal computer again] 

: s (E-v"kb) [look for pair with 2nd > 16k + # of bits read] 

/s +s [Found it! Output 1st and halt] 

/="!+v (Lt''b*@()) [Read another bit from program tape] 

/="?+v (L*ltb) [Increase depth/time limit] 

"? [Surprise, formal system halts, so we do too] 

(L()()) [Initially, depth limit and no bits read] 

[ 

,#,','((xy)(ll)) 
] 

,#,','(x()) 



godelS.lisp 



[[[ 

Show that a formal system of complexity N 
can't determine more than N + c bits of Omega, 
and also determine c. 

Formal system is a never halting lisp expression 
that outputs lists of the form (lOXOXXXXlO) . 
This stands for the fractional part of Omega, 
and means that these 0,1 bits of Omega are known. 
X stands for an unknown bit. 
]]] 

[Count number of bits in an omega that are determined.] 
& (Cw) /.w() - /=0+w'(l) /=l+w'(l) 
(C-w) 

(C (XXX)) 
(C'(IXX)) 
(C'(IXO)) 
(C'(llO)) 

[Merge bits of data into unknown bits of an omega.] 
& (Mw) /.w() * /=0+wO /=l+wl @ 
(M-w) 

[Test it.] 

++?0 ' : (Mw)/.w()*/=0+wO/=l+wl@(M-w) (M' (OOXOOXOOX) ) ^111) 
++?0 ' : (Mw)/.w()*/=0+wO/=l+wl@(M-w) (M' (11X11X111)) '(00) 

[(<xy) tells if x is less than y.] 
& (<xy) /.X /.yOl 

/.yO«-x-y) 
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(<'(11)'(11)) 
(<' (11) ' (111)) 
«'(111)'(11)) 

[ 

Examine omegas in list w to see if in any one of them 
the niimber of bits that are determined is greater than n. 
Returns to indicate not found, or what it foimd. 

] 

& (Ewn) /.w /(<n(C+w)) +w (E-wn) 

(E'((00)(000))'()) 

(E'((00)(000))'(1)) 

(E' ((00) (000)) '(ID) 

(E' ((00) (000)) ' (111)) 

(E' ((00) (000)) '(1111)) 

++?0' !% 

[Count number of bits in an omega that are determined.] 
: (Cw) /.w() - /=0+w'(l) /=l+w'(l) 
(C-w) 

[Merge bits of data into unknown bits of an omega.] 
: (Mw) /.w() * /=0+wO /=l+wl @ 
(M-w) 

[(<xy) tells if x is less than y.] 
: «xy) /.X /.yOl 

/.yO«-x-y) 
[Over-write real definition for test.] 
: «xy) 1 

[ 

Examine omegas in list w to see if in any one of them 
the number of bits that are determined is greater than n. 
Returns to indicate not found, or what it foimd. 

] 

: (Ewn) /.w /(<n(C+w)) +w (E-wn) 



godelS.lisp 



51 



[Parameter in proof] 
k -'(11111) 
k ~ 'kk 
k ~ 'kk 
k ~ 'kk 
k ~ 'kk 

[Main Loop: t is depth limit (time), b is bits read so far.] 
: (Ltb) 

: V ~?~t'!%~b [Run universal computer again.] 
: s (E-v'kb) [Look for an omega with > 

(16k + # of bits read) bits determined.] 
/s (Ms) [Found it! Merge in undetermined bits, 

output result, and halt.] 
/="!+v (Lt'b*@()) [Read another bit from program tape.] 
/="?+v (L*ltb) [Increase depth/time limit.] 

"? [Surprise, formal system halts, 

so we do too.] 

(LOO) [Initially, depth limit 

and no bits read.] 

,'(1X0) [Toy formal system with only one theorem.] 



(0) [Missing bit of omega that is needed.] 
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Part III 
The Software 
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/* lisp.c: self-contained high-speed LISP interpreter */ 
/* 

The storage required by this interpreter is 8 bytes times 
the symbolic constant SIZE, which is 8 * 10,000,000 = 
80 megabytes. To run this interpreter in small machines, 
reduce the #define SIZE 10000000 below. 

To compile, type 

cc -0 -olisp lisp.c 
To run interactively, type 

lisp 

To run with output on screen, type 

lisp <test.lisp 
To run with output in file, type 

lisp <test.lisp >test.run 

Reference: Kernighan & Ritchie, 

The C Programming Language, Second Edition, 

Prentice-Hall, 1988. 

*/ 

#include <stdio.h> 
#include <time.h> 

#define SIZE 10000000 /* numbers of nodes of tree storage */ 
#define LAST_ATOM 128 /* highest integer value of character */ 
#define nil 128 /* null pointer in tree storage */ 
#define question -1 /* error pointer in tree storage */ 
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#define exclamation -2 /* error pointer in tree storage */ 
#define infinity 999999999 /* "infinite" depth limit */ 

long hd[SIZE+l] , tl[SIZE+l]; /* tree storage */ 

long next = nil; /* list of free nodes */ 

long low = LAST_AT0M+1; /* first never-used node */ 

long vlst [LAST_AT0M+1] ; /* bindings of each atom */ 

long tape; /* Turing machine tapes */ 

long display; /* display indicators */ 

long outputs; /* output stacks */ 

long q; /* for converting expressions to binary */ 

long col; /* column in each 50 character chunk of output 

(proceeded by 12 char prefix) */ 
time_t timel; /* clock at start of execution */ 
time_t time2; /* clock at end of execution */ 

long ev(long e) ; /* initialize and evaluate expression */ 

void initialize_atoms(void) ; /* initialize atoms */ 

void clean_env(void) ; /* clean environment */ 

void restore_env(void) ; /* restore dirty environment */ 

long evaldong e, long d) ; /* evaluate expression */ 

/* evaluate list of expressions */ 

long evalstdong e, long d) ; 

/* bind values of arguments to formal parameters */ 

void binddong vars, long args) ; 

long at (long x) ; /* atomic predicate */ 

long jndong x, long y) ; /* join head to tail */ 

long pop (long x) ; /* return tl & free node */ 

void frdong x) ; /* free list of nodes */ 

long eqdong x, long y) ; /* equal predicate */ 

long cardinality (long x) ; /* number of elements in list */ 

long appenddong x, long y) ; /* append two lists */ 

/* read one square of Turing machine tape */ 

long getbit (void) ; 

/* read one character from Turing machine tape */ 
long getchr (void) ; 

/* read expression from Turing machine tape */ 
long getexpdong top) ; 

void putchrdong x) ; /* convert character to binary */ 
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void putexpClong x) ; /* convert expression to binary */ 
long outCchar *x, long y) ; /* output expression */ 
void out2(long x) ; /* really output expression */ 
void out3(long x) ; /* really really output expression */ 
long chr2(void); /* read character - skip blanks, 

tabs and new line characters */ 
long chr(void); /* read character - skip comments */ 
long indong mexp, long rparenokay) ; /* input s-exp */ 

mainO /* lisp main program */ 
{ 

char name_colon[] = "X:"; /* for printing name: def pairs 

timel = time(NULL); /* start timer */ 
printf ("lisp.c\n\nLISP Interpreter Run\n"); 
initialize_atoms() ; 

while (1) { 

long e , f , name , def ; 
printf ("\n") ; 

/* read lisp meta-expression, ) not okay */ 
e = in(l,0); 

/* flush rest of input line */ 
while (putchar (getchar ) != '\n'); 
printf ("\n") ; 
f = hd[e] ; 
name = hd [tl [e] ] ; 
def = hd[tl[tl[e]]] ; 
if (f == { 
/* definition */ 
if (at (name)) { 

/* variable definition, e.g., & x ' (abc) */ 

def = out (" expression" , def ) ; 

def = ev(def ) ; 
} /* end of variable definition */ 
else { 

/* function definition, e.g., & (Fxy) *x*y() */ 
long var_list = tl [name] ; 
name = hd[name] ; 
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def = jn('&' , jn(var_list, jn(def ,nil))) ; 

} /* end of function definition */ 

naiiie_colon[0] = name; 

out(name_colon,def ) ; 

/* new binding replaces old */ 

vl St [name] = jn(def,nil); 

continue ; 
} /* end of definition */ 
/* write corresponding s-expression */ 
e = out (" expression" ,e) ; 
/* evaluate expression */ 
e = out ( "value " ,ev(e)) ; 

} 

} 



long evdong e) /* initialize and evaluate expression */ 

{ 

long d = infinity; /* "infinite" depth limit */ 
long v; 

tape = jn(nil,nil); 

display = jn('Y',nil); 

outputs = jn(nil,nil); 

V = eval(e,d) ; 

if (v == question) v = '?' ; 

if (v == exclamation) v = ' ! ' ; 

return v; 

} 

void initialize_atoms (void) /* initialize atoms */ 
{ 

long i; 

for (i = 0; i <= LAST_ATOM; ++i) { 

hd[i] = tl [i] = i; /* so that hd & tl of atom = atom */ 
/* initially each atom evaluates to self */ 
vlst[i] = jn(i,nil); 
} 

} 



long jndong x, long y) /* join two lists */ 
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{ 

long z; 

/* if y is not a list, then jn is x */ 
if ( y != nil && at(y) ) return x; 

if (next == nil) { 
if (low > SIZE) { 
printf ("Storage overf low! \n") ; 
exit(O) ; 
} 

next = I0W++; 
tl[next] = nil; 
} 

z = next; 
next = tl [next] ; 
hd[z] = x; 
tl[z] = y; 

return z; 

} 

long pop (long x) /* return tl & free node 
{ 

long y; 
y = tl [x] ; 
tl [x] = next ; 
next = x; 
return y; 

} 

void frdong x) /* free list of nodes */ 
{ 

while (x != nil) x = pop(x); 

} 

long at (long x) /* atom predicate */ 
{ 

return ( x <= LAST.ATOM ) ; 
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y 

long eqClong x, long y) /* equal predicate */ 
{ 

if (x == y) return 1 ; 
if (at(x)) return 0; 
if (at(y)) return 0; 

if (eq(hd[x] ,hd[y])) return eq(tl [x] ,tl [y] ) ; 
return 0; 

} 

long evaldong e, long d) /* evaluate expression */ 
{ 

/* 

e is expression to be evaluated 

d is permitted depth - integer, not pointer to tree storage 
*/ 

long f, V, args, x, y, z, vars, body; 

/* find current binding of atomic expression */ 
if (at(e)) return hd[vlst[e]]; 

f = eval(hd[e] ,d) ; /* evaluate function */ 

e = tl[e]; /* remove f miction from list of arguments */ 

if (f < 0) return f; /* function = error value? */ 

if (f == '\") return hd[e] ; /* quote */ 

if (f =='/'){/* if then else */ 
V = eval(hd[e] ,d) ; 
e = tl[e] ; 

if (v < 0) return v; /* error? */ 
if (v == '0') e = tl[e] ; 
return eval(hd[e] ,d) ; 
} 

args = evalst(e,d); /* evaluate list of arguments */ 
if (args < 0) return args; /* error? */ 
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X = hd[args] ; /* pick up first argument */ 

y = hd[tl[args]] ; /* pick up second argument */ 

z = hd[tl[tl[args]]] ; /* pick up third argument */ 



switch 
case 
case 
case 



(f) { 

{fr(args); return getbitO;} 
{fr(args); return getexpC 'Y' ); } 
{fr(args) ; 

V = q = jn(nil,nil); putexp(x); return popCv);} 
{fr(args) ; return hd[x] ;} 
{fr(args); return tl[x];} 
{fr(args); return (at(x) ? '1' : '0');} 
{fr(args) ; hd [outputs] = jn(x,hd [outputs] ) ; 
return (hd[display] == 'Y' ? out ( "display" ,x) : x);} 
{fr(args) ; return out ( "display " ,x) ;]• 
{fr(args); return (eq(x,y) ? '1' : '0');} 
{fr(args); return jn(x,y);} 
{fr(args) ; 

return append ( (at (x)?nil:x) , (at(y)?nil:y)) ;} 



if (d == 0) {fr(args); return question;} /* depth exceeded 

-> error! */ 

d — ; /* decrement depth */ 

if (f == ' ! ') { 
fr(args) ; 

clean_env() ; /* clean environment */ 

V = eval(x,d) ; 

restore_env() ; /* restore unclean environment */ 

return v; 

} 

if (f == '?') { 
fr(args) ; 

X = cardinality (x) ; /* convert s-exp into number */ 

clean_env() ; 

tape = jn(z,tape); 

display = jn( 'N' , display) ; 
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outputs = jn(nil , outputs) ; 

V = eval(y,(d <= x ? d : x)); 
restore_env() ; 

z = hd [outputs] ; 
tape = pop(tape); 
display = pop (display) ; 
outputs = pop(outputs) ; 

if (v == question) return (d <= x ? question : jn('?',z)); 
if (v == exclamation) return jn('!',z); 
return jn(jn(v,nil) ,z) ; 
} 

f = tl[f] ; 
vars = hd [f ] ; 
f = tl[f] ; 
body = hd[f] ; 

bind(vars,args) ; 
fr(args) ; 

V = eval(body,d) ; 

/* unbind */ 

while (! at (vars)) { 

if (at(hd[vars] ) ) 

vlst [hd [vars] ] = pop (vlst [hd [vars]] ) ; 

vars = tl [vars] ; 

} 

return v; 

} 

void clean_env(void) /* clean environment */ 
{ 

long i; 

for (i = 0; i <= LAST_ATOM; ++i) 

vlst [i] = jn(i , vlst [i] ) ; /* clean environment */ 

} 
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void restore_env(void) /* restore unclean environment */ 
{ 

long i; 

for (i = 0; i <= LAST_ATOM; ++i) 

vlst [i] = popCvlst [i] ) ; /* restore unclean environment */ 

} 

long cardinality (long x) /* number of elements in list */ 
{ 

if (at(x)) return (x == nil ? : infinity); 
return l+cardinality(tl [x] ) ; 

} 

/* bind values of arguments to formal parameters */ 

void binddong vars, long args) 

{ 

if (at (vars)) return; 
bind(tl [vars] ,tl [args] ) ; 
if (at (hd [vars])) 

vlst [hd [vars] ] = j n(hd [args] , vlst [hd [vars]] ) ; 

} 

long evalstdong e, long d) /* evaluate list of expressions 
{ 

long X, y; 

if (at(e)) return nil; 

X = eval (hd [e] , d) ; 

if (x < 0) return x; /* error? */ 

y = evalst(tl[e] ,d) ; 

if (y < 0) return y; /* error? */ 

return jn(x,y) ; 

} 

long append (long x, long y) /* append two lists */ 
{ 

if (at(x)) return y; 

return jn(hd[x] , append (tl [x] ,y)) ; 

} 
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/* read one square of Turing machine tape */ 

long getbit (void) 

{ 

long x; 

if (at (hd [tape] ) ) return exclamation; /* tape finished ! */ 
X = hd[hd[tape]] ; 
hd[tape] = tl[hd[tape]] ; 
return (x == '0' ? '0' : '10; 

} 

/* read one character from Turing machine tape */ 

long get chr (void) 

{ 

long c, b, i; 
c = 0; 

for (i = 0; i < 7; ++i) { 
b = getbit ; 

if (b < 0) return b; /* error? */ 

c=c+c+b- '0'; 

} 

/* nonprintable ASCII -> ? */ 

return (c > 31 && c < 127 ? c : '?'); 

} 

/* read expression from Turing machine tape */ 

long getexpdong top) 

{ 

long c = getchrO, first, last, next; 
if (c < 0) return c; /* error? */ 

if (top == 'Y' && c == ')') return nil; /* top level only */ 
if (c != '(') return c; 
/* list */ 

first = last = jn(nil,nil); 

while ((next = getexp('NO) != ')') { 

if ( next < ) return next; /* error? */ 

last = tl[last] = jn(next ,nil) ; 

} 

return pop(first); 

} 
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void putchrdong x) 
{ 

q = tl[q] = jn(( x 
q = tl[q] = jn(( x 
q = tl[q] = jn(( X 
q = tl[q] = jn(( X 
q = tl[q] = jn(( x 
q = tl[q] = jn(( x 
q = tl[q] = jn(( x 

} 

void putexpClong x) 
{ 

if ( at(x) && x != 
putchrC (') ; 

while (!at(x)) { 
putexp (hd [x] ) ; 
X = tl [x] ; 
} 

putchrC) ') ; 

} 

long out (char *x, long y) /* output expression */ 
{ 

printf ("y.-12s",x); 

col =0; /* so can insert \n and 12 blanks 

every 50 characters of output */ 

out2(y) ; 
printf ("\n") ; 
return y; 

} 

void out2(long x) /* really output expression */ 
{ 

if ( at(x) && X != nil ) {out3(x); return;} 
outSCC); 



/* convert character to binary */ 
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/* convert expression to binary */ 
nil ) {putchr(x); return;} 
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while (!at(x)) { 
out2(hd [x] ) ; 
X = tl [x] ; 
> 

outSC)'); 



void outSdong x) /* really really output expression */ 
{ 

if (col++ == 50) {printf ("\n%-12s"," "); col = 1;} 
putchar (x) ; 

} 

long chr2(void) /* read character - skip blanks, 

tabs and new line characters */ 

{ 

long c; 
do { 

c = getchar ; 
if (c == EOF) { 

time2 = time (NULL) ; 
printf ( 

"End of LISP Rim\n\nElapsed time is %.0f seconds. \n", 
dif f time (time2 , timel) 

); 

exit(O); /* terminate execution */ 
} 

putchar (c) ; 

> 

/* keep only non-blank printable ASCII codes */ 
while (c >= 127 I I c <= 32) ; 
return c; 

} 

long chr(void) /* read character - skip comments */ 
{ 

long c; 
while (1) { 
c = chr2() ; 
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if (c != 'L') return c; 

while (chr() != ; /* comments may be nested */ 

} 

} 

long indong mexp, long rparenokay) /* input m-exp */ 
{ 

long c = chrO ; 

if (c == ')') if (rparenokay) return ')'; else return nil; 
if (c == '(') { /* explicit list */ 

long first, last, next; 

first = last = jn(nil,nil); 

while ((next = in(mexp,l)) != ')') 

last = tl[last] = jn(next ,nil) ; 

return pop(first); 

} 

if (!mexp) return c; /* atom */ 
if (c == "") return in(0,0); /* s-exp */ 
if (c == ':') { /* expand "let" */ 
long name , def , body ; 
name = in(l ,0) ; 
def = in(l,0); 
body = in(l,0) ; 
if (! at (name)) { 
long var_list; 
var_list = tl [name] ; 
name = hd[name] ; 
def = 

jn('\' ' , jn(jn('&' , jn(var_list, jn(def ,nil))) ,nil)) ; 
} 

return 
jn( 

jn('\' ' , jn(jn('&' , jn(jn(name,nil) , jn(body,nil) ) ) ,nil)) 
jn(def ,nil) 

); 
} 

switch (c) { 

case ' ' : case '°/o' : 
return jn(c,nil); 
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case '+' : case case ' . ' : case '\'': 

case ',': case case '#': case 

return jn(c, jn(in(l,0) ,nil)) ; 
case '*': case '=' : case 'k': case 

return jn(c, jn(in(l,0) , jn(in(l,0) ,nil))) ; 
case '/': case ' : ' : case : 

return jn(c, jn(in(l,0) , jn(in(l,0) , jn(in(l,0) ,nil)))) ; 
default : 

return c; 

} 



